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ABSTRACT 

Due to their low mass and luminosity, M dwarfs are ideal targets if one hopes to find low-mass planets similar to Earth by using 
the radial velocity (RV) method. However, stellar magnetic cycles could add noise or even mimic the RV signal of a long-period 
companion. Following our previous work that studied the correlation between activity cycles and long-term RV variations for K 
dwarfs we now expand that research to the lower-end of the main sequence. Our objective is to detect any correlations between long- 
term activity variations and the observed RV of a sample of M dwarfs. We used a sample of 27 M-dwarf s with a median observational 
timespan of 5.9 years. The cross-correlation function (CCF) with its parameters RV, bisector inverse slope (BIS), full-width-at-half- 
maximum (FWHM) and contrast have been computed from the HARPS spectrum. The activity index have been derived using the Na i 
D doublet. These parameters were compared with the activity level of the stars to search for correlations. We detected RV variations 
up to ~5 m s~^ that we can attribute to activity cycle eff'ects. However, only 36% of the stars with long-term activity variability appear 
to have their RV aff'ected by magnetic cycles, on the typical timescale of ~6 years. Therefore, we suggest a careful analysis of activity 
data when searching for extrasolar planets using long-timespan RV data. 
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1. Introduction 

The majority of extrasolar planets discovered so far were either 
detected or confirmed via the radial velocity methocfl. This tech- 
nique measures the Doppler efiTect due to the wobble of the star 
around the centre of mass of the star-planet system. As an indi- 
rect method, it is sensitive to stellar sources of noise like oscilla- 
tions, granulation, rotating active regions, and magnetic cycles. 

Stellar oscillations and granulation induce radial veloc- 
ity variations on timescales up to some hours and can eas- 
ily be suppressed if a op timized observational strategy is used 
(iDumusque et al.|[201 Ibh . However, this is not true for the case 
of rotationally modulated active regions which have longer 
timescales. The eflTect of rotating active regions on line profiles, 
and hence RV measurement s, can easily hide or mimic the signal 
of orbiting companions (e .g. lSaar & Donahue! 1997l:ISantos et al.l 
120001: lOueloz et al. 11200 Ih . In some cases, these eflTects can be di- 
agnosed e.g. by the anti-correlation between instan taneous mea- 
surements of RV and the bisector inverse slope (lOueloz et al.l 
[2OQI : Boisse et al. 2011) and corrected by su btracting the anti- 
correlation slope t o the RV measurements (e.g. lMelo et al .1120071: 
iBoisse et"aDl2009D or by fitting an extra Keplerian orbit to the 



* Based on observations made with the HARPS instrument on the 
ESO 3.6-m telescope at La Silla Observatory under programme ID 
072.C-0488(E2 

^ cf . jhttp://exoplanet.eu71 



RV data with the pe riod detected in the act ivity time- series (e.g . 
iBonfils et al.l I20Q7I: iForveille et al.l l2009h . lOueloz et al.l (l2009l) 
and lBoisse et al.l ( 2011 ) proposed to fit the activity signal with 
sinusoids of different periods: the rotation one and its harmon- 
ics. Saar & Fischer ( 2000) used a diflTerent technique to correct 
for long-term activity induced RV. These authors used the slope 
of the S iR-RV correlation to remove the activity influence from 
the radial velocity signal. 

Long-term stellar magnetic cy cles can also be a so urce of 
noise for precise RV measurements. iKiirster et al.l (l2003l) studied 
the correlation between radial velocity and the Ha index for the 
nearby Barnard's star (G1699). They found an anti-correlation 
between velocity and the activity index with a correlation co- 
eflftcient of p = -0.50. The authors concluded that activity as 
measured by the hydrogen line was producing a blueshift of the 
photospheric absorption lines. 

By using simulations, iMeunier et al.l (1201 Ol) showed that 
magnetic activity cycles can induce radial-velocity variations for 
the case of the Sun as seen edge-on with amplitudes that can 
reach the ~10 m s"^ level. These kind of variations in a star with 
a periodic activity cycle could be able to mimic the signal of a 
long-period extra- solar planet. 

In this context. ISantos et aP (1201 Oh used a sample of 8 solar- 
type stars with the aim of studying if a correlation between long 
term activity and RV variations exists. The long-term variations 
were detected on the S mw, Ha, and Hei indices and the BIS, 
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FWHM and contrast CCF parameters but only two stars showed 
correlations with RV stronger than p > |0.75|: one positive corre- 
lation and one anti-correlation. The authors concluded then that 
the possible amplitudes of induced RV variations for the early- 
K dwarfs was low, at the ~1 ms"^ level similar to the HARPS 
precision. 

Using a larger sample of ar ound 300 sta r s from the HARPS 
FGK high precision program, iLovis et aL 1 (l20TTh studied the 
correlation between long-term activity varia t ions w ith radial- 
velocity in a similar fashion as ISantos et"an (1201 Ol) did. They 
found a correlation between the slope of the RV-activity index 
(/?^^) correlation, effective temperature (Tq^), and metallicity 
([Fe/H]). The slope is weaker for late-type dwarfs than for early 
ones and therefore, it seems that the RVs of later-K dwarfs are 
less affected by magnetic cycles than the RVs of early-G dwarfs. 
Therefore, the influence of long-term activity could be corrected 
if the activity level, effective temperature and metallicity of the 
star are known by using the slope of the RV-log(7?^^) relation. 

Hints of long t erm RV varia t ions p roduced by activity cy- 
cles were found bv lMoutou et al.l (l201lh in the stars BD-1 14672 
and HIP21934, with pe riods of 1692 and 1 100 days respectively. 
ISegransan et al.l (l201ll) also found a long term period of -500 
days on the radial velocity of HD 104067 (K2V). The activity in- 
dex of this star, hosting a 5 5 -day period Neptune-like planet, 
was found to have a correlation with the RV residuals (with 
o-Rv = 4.6 m s"^) after the planetary signal was removed. 

More recently, three exoplanets were discovered in three 
early-K stars with magnetic cycl es by correcting the ac tivity 
signals in the radial velocity data (iDumusque et al.|[2bl lah . The 
planets were found by fitting simultaneously two keplerians: one 
for the planet and one for the magnetic cycle. All the parame- 
ters of the keplerian fitting the cycle, except the amplitude, were 
fixed to the ones obtained when fitting the activity index only. 
This proves that (i) magnetic activity cycles of stars can influ- 
ence radial velocity and hide the signal of long-period planets, 
and (ii) a correction of the long-term activity-induced RV is pos- 
sible and can be used to recover the embedded signal of a planet. 

In iGomes da Silva et aP (l201ll hereafter Paper I) we com- 
pared the long-term activity variations using four activity indices 
for a sample of M-dwarf stars from the HARPS program. We 
arrived to the conclusion that the Na i index was the most appro- 
priate of the four to study the activity of these type of stars. We 
will now use in this paper the Na i index as the proxy of activity 
with the aim of comparing long-term activity to RV and some 
other relevant CCF parameters using the same san iple of stars. 
This p aper will therefore extend the first study in ISantos et"aD 
(I2OIOI) to the case of early-M dwarfs. 

This paper is organised as follows: in Sect. [2l we present our 
sample and observation log and explain our data analysis, in Sect 
[3] a selection of a subsample of stars with long-term activity vari- 
ability is done, in Sect. |4] we determine which stars have hints of 
periodic magnetic cycles, in Sect. \5\ we compare long-term ac- 
tivity with radial velocity and the CCF parameters, in Sect. 
we examine individual stars with strong activity-RV correlation 
and other interesting cases, and finally in Sect. [7] we draw our 
conclusions from the present work. 

2. Sample and observations 

The sample we are using comes from the HARPS M-dwarf 
pla net search program that started in 2003 and ended in 2009 
(see lBonfils et aLluOllI) . We used this sample in Paper I to study 
four known chromospheric activity indices and to select stars 
with long-term activity variability. However, we now use more 



data that were taken in 2010 as part of the lBonfils et al.l (l201ll) 
program extension. Therefore, We decided to redo the same 
analysis as in Paper I using the new data to detect any new cases 
of activity variability that could arise from more data points. 

We obtained simultaneous radial- velocity, BIS, FWHM, 
contrast, and the Nai activity indicator. The median radial- 
velocity error of the nightly averaged data was 1.2 ms"^. 
HARPS is capab le of more precise measurements (e.g. 
iMayor et"an l201 ll) but our sample includes dim stars which 
make it more difficult to get high signal-to-noise ratios (SNR) 
and therefore higher-precision RVs. To measure activit}Q we 
used the index based on the Na i D 1 and D 2 lines in a same 
way as in Paper I (see also lDiaz et al The CCF parame- 

ters were used because of their potential as co mplementary long- 
term activity proxies (e.g. lSantos et al.|[2bl0l) . 

Because of an instrumental drift detected on HARPS that 
affected the FWHM and contrast of the CCF, these parameters 
were corrected using the following expressions: 

FWHM^orr = FWHM -h 5.66 • 10"^Z) - 1.77 • 10"^Z)^ (1) 
contrastcorr = contrast -h 4.59 • lO'^D - 2.75 • 10"^D^ (2) 

where D is (BJD - 2454500) days. However, this correction is 
very small and is a long-term drift. For example the FWHM drift 
is only 0.1% in five years. 

Since this sample includes some stars very close to the Sun, 
these stars may have significant secular acceleration that could 
produce a trend on the observed radial- velocity. We corrected all 
stars for this effect using the proper motio ns and para llaxes that 
were retrieved from the Hipparcos catalog (lES All 19971) . The sec- 
ular accelerations were calculated following the description in 
IZechmeister et al.l (l2009l) and the radi al-velocities subsequently 
corrected (see also lBonffis et al.|[201ll) . 

Some of the stars in this sample are planet hosts, namely 
Gl 176, Gl 433, Gl 436, Gl 581, Gl 667 C, G1674 G183 2, Gl 849, 
Gl 876, HIP12961, and HIP85647 (iBonffis et al.ll201lh . Because 
we are only interested in variations due to activity, these plane- 
tary signals were subtracted from the observed radial-velocities 
by fitting keplerian functions based on the published orbital pa- 
rameters. 

The selection criteria was the same as in our previous pa- 
per. We nightly averaged and binned the data into 150-day bins 
to average out short time- scale variability, since we are only 
concerned with long-term variations. Only bins with more than 
three nights and stars with at least four bins were selected. The 
errors used for each bin are the statistical errors on the aver- 
age, (Ti = cr/ Va^, where cr is the root-mean- square (rms) of the 
nightly averaged data in each bin, and N is the number of nights 
included in the bin. 

We also decided to select only stars with at least three years 
of observations, which resulted in GJ 361, GJ 2049, and GJ 3218 
being discarded from the sample. This resulted in a sample of 
27 stars which passed the selection criteria, whose basic param- 
eters are presented in Table [T] These stars have spectral types 
in the range M0-M5.5 V and have V magnitudes between 7.34 
and 1 1.81. In the same table we also present the observation log, 
with information about the number of nights of observation, the 
timespan in days and years, the average signal-to-noise ratio at 
spectral order 56 (the order of the Nai doublet, -5893 A) and the 
average radial- velocity error. The timespan of the observations 
range from 3.3 to 7.2 years (with a median value of 5.9 years). 

^ From now on we will use the terms activity, Na i, and sodium lines 
interchangeably. 
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Table 1: Basic parameters and observation log of the sample. 



Star Other names Sp. Type^ V-T BJD.tan 2400000 # nights < S IN > < (T,(K) > 







[mag] 


[mag] 


[days] 




[days] 


(order 56) 


[ms-^] 


Gil 


M3V 


8.57 


2.13 


52985.60 


45 


2063 


95.6 


0.7 


Gl 176^ 


M2.5V 


9.97 


2.25 


52986.71 


71 


2146 


47.6 


1.2 


G1205 


M1.5V 


7.92 


2.08 


52986.73 


75 


1400 


108.5 


0.8 


G1273 


M3.5V 


9.89 


2.71 


52986.77 


62 


2140 


57.3 


0.8 


G1382 


M2V 


9.26 


2.18 


52986.84 


30 


1188 


66.7 


1.1 


G1393 


M2V 


9.76 


2.26 


52986.86 


29 


1897 


67.0 


1.0 


G1433^ 


M2V 


9.79 


2.15 


52989.84 


61 


2245 


56.6 


1.2 


G1436^ 


M2.5V^ 


10.68 


2.02 


53760.83 


115 


1532 


32.9 


1.4 


G1479 


M3V 


10.64 


1.90 


53158.55 


58 


1413 


43.3 


1.2 


G1526 


M1.5V 


8.46 


2.07 


53158.60 


34 


1843 


101.6 


0.8 


G1551 


Prox Cent M5.5V 


11.05 


3.62 


53152.60 


42 


2140 


21.1 


1.3 


G1581^ 


M2.5V 


10.57 


2.53 


53152.71 


194 


2306 


36.4 


1.2 


G1588 


M2.5V 


9.31 


2.40 


53152.75 


32 


2200 


75.4 


0.8 


G1667C^ 


M2V 


10.22 


2.08 


53158.76 


173 


2201 


42.9 


1.3 


G1674^ 


M3V 


9.36 


2.40 


53158.75 


44 


1574 


70.6 


0.8 


G1680 


M1.5V 


10.14 


2.27 


53159.71 


35 


2269 


41.2 


1.4 


Gl 699 


Barnard's star M4V 


9.54 


2.52 


54194.89 


32 


1291 


34.3 


0.8 


G1832^ 


MIV 


8.67 


2.18 


52985.52 


59 


2424 


88.9 


0.8 


G1849^ 


M3V 


10.42 


2.50 


52990.54 


49 


2423 


41.9 


1.3 


G1876^ 


M3.5V 


10.17 


2.77 


52987.58 


67 


2508 


42.6 


1.1 


G1877 


M2.5V 


10.55 


2.43 


52857.83 


40 


2636 


36.5 


1.2 


G1887 


M2V 


7.34 


2.02 


52985.57 


63 


1407 


131.3 


0.8 


G1908 


MIV 


8.98 


2.04 


52986.58 


66 


2511 


76.0 


0.9 


HIP 12961^ 


MOV" 


9.7 


1.64 


52991.63 


46 


2226 


46.5 


3.0 


HIP 19394 


M3.5V^ 


11.81 


2.5 


52942.80 


35 


2495 


23.3 


2.0 


HIP38594 


M^ 


9.96 


1.64 


52989.79 


17 


2229 


41.1 


2.0 


HIP85647^ 


GJ676A MOV^ 


9.59 


1.85 


53917.75 


38 


1520 


36.6 


2.4 



Notes. The average radial-velocity errors, < cr,(y^) >, are taken after removal of the planetary companion's signals, (t) Stars with published 
planetary companions. 

References. (a) Bon fils et al] (l201 1 | ) unless individually specified; (b) iHawlev et all (119961) : (c) iKoen et all (l201Qh : (d) Simbad 
( |http :// simbad. u- strasbg . fr/simbad/ ) ; (e ) ]ESAl ( [T997h 



The average errors on radial velocity per star for the nightly av- 
eraged data vary between 0.8 and 3.0 m s"^ 

3. Nai index variability 

In Paper I we compared four activity indices, namely 5'caii, Ha, 
Na I, and He l, and arrived to the conclusion that the Na I in- 
dex is the best to follow the long-term activity of IM-dwarf stars. 
We also tested our sample for variability and found that Gl 1 , 
Gl 273, Gl 433, Gl 436, Gl 581, Gl 588, Gl 667 C, Gl 832, Gl 849, 
Gl 877, Gl 908, and HIP85647 showed significant {P{F) < 0.05) 
long-term activity variability on at least two indices. However, 
we are now using the latest data with new measurements from 
2010 (previously we used data covering the years 2003-2009) 
and consequently other stars might now show statistical variabil- 
ity on long time- scale activity that was not detected before. We 
therefore decided to repeat the variability F-tests, this time using 
only the Na I index. 

The Na i activity proxy was determined as explained in Paper 
I, by measuring the flux in the centre of the sodium Dl and D2 
lines in comparison to the flux on two reference bands. This in- 
dex is not calibrated to the bolometric flux of the stars and there- 
fore it will depend on the eflTective temperature of each star (see 
discussion on Sect. 6.6 of Paper I). As a consequence, this in- 
dex cannot be used to compare the activity of stars with diflTerent 
eflTective temperatures. Nevertheless, it can be used to detect ac- 
tivity variability with time for a given star. 



In a similar fashion to what was done in Paper I, we asked 
whether the long time-scale variations observed in our data were 
of statistical significance. As previously done in Paper I, we 
used the F-test with an F- value as F = cr^/ < ct/ >^ where 
(Te is the standard deviation and < ct/ > the average of the 
error on the me a n of the binned parame t er for each star (se e 
lEndl et al.1 [20021: IZechmeister et al.l l2009l: iBonfils et al.l l201 ih . 
This F-test will give the probability that the observed standard 
deviation can be explained by the random scatter due to the in- 
ternal errors. Therefore, a low value of P{F) will discriminate 
the stars which have significant variability not justified by the 
internal errors. The results of these tests are shown in Table O 
Bold values indicate probabilities lower than 0.05 (95% signifi- 
cance level). The F- value used to calculate the probabilities was 
F = (Te(Na l)^/ < CT/ (Na l) >^ where (Te stands for the rms of the 
(binned) data and < ct/ > is the average of the error on the mean. 
Also shown is the number of bins Nbins, indication if the star 
passed the variability tests in Paper I, SeL, the radial-velocity's 
rms < (TeiVr) > and average error on the mean < cri(Vr) >, and 
the average value of the Nai index, < Na I >. 

Fourteen out of the 27 stars show long-term activity vari- 
ability with P(F) < 0.05, representing 52% of the sample. 
These stars are G1273, G1433, G1436, G1526, G1581, G1588, 
G1667C, G1699, G1832, G1849, G1876, G1877, G1908, and 
HIP85647. These are the same stars that passed the tests in Paper 
I plus Gl 526, Gl 699, and Gl 876, and less Gl 1. Four more stars, 
G1680, G1887, HIP19394, and HIP38594, have 0.05 < P(F) < 
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0.1, which if they are considered increases the percentage of 
stars with variability to 67%. In Paper I we found 17 stars out 
of 30 showing variability in the Nai index with P(F) < 0.1, 
representing 57 % of the sample. This represents an increase of 
variable stars by 10% from our previous work just by adding one 
more year of observations, and thus, more of these stars could 
show long-term variability if a longer time-span is considered. 
We should note that Gl 1 passed the variability test on Paper I 
but not here. This is because on Paper I this star passed the test 
due to its variability in the 5" can and Ha indices, while not show- 
ing significant variability in Na i. Since we used only this index 
now and Gl 1 have not had more data points in 2010, this star 
now failed to pass the F-test based just on Nai. 

By using a sample of old FGK stars in the s olar neighbour- 
hood studied for up to 7 years jLovis et al.l (1201 ll) found that 6 1 % 
have a magnetic activity cycle. Although we cannot confirm the 
cyclic nature of our detected long- term activity varia tions, our 
results are compatible with those of iLovis et al.l (1201 

From now on we will consider only these 14 stars for the rest 
of the study, unless specified. 

4. Activity cycle fits 

Some stars show clear long-term variation of their activity, well 
visible with a general trend (increase or decrease of activity 
during several observational seasons). These include stars that 
have not passed the variability F-tests of Sect. [3l There is evi- 
dence that the periodic activity cycles of M-dwarf st ars can be 
well fitted by a sinusoidal function (see .Cincunegui et al. 20071; 
iDiaz et al.1 r2007bl: iBuccino et al.l l201 ih . Given the usually re- 
duced number of points in our study, fitting other functions (e.g. 
Keplerian) would not bring further constraints since it would in- 
crease the number of free fitting parameters close to the number 
of data points. We therefore tried to fit a sinusoidal signal to Na i 
data to obtain the periods of such cycles (minimum periods for 
the cases of data covering less than one period). The fitted sinu- 
soidal signal was of the form 



Na isine = A sin 



27r 

P cycle 



(3) 



where A is the semi- amplitude in units of Na i, P cycle the activity 
cycle period in days, t the observation epoch in days, the phase, 
and 6 a constant off'set. The quality of the fit was obtained by 
a p-value, which gives the probability that the data is better fit 
by a straight line than by a sinusoidal signal. This p-value was 
calculated via an F- value given by 



sine 



(4) 



where line and sine identify if the model used was the linear fit 
or sinusoidal fit, respectively, and N is the number of data points. 
The chi- squared function used was 



X 



Nai„ - Nai„,, 



cr/(Nai)„ 



, for 1 < < AA, 



(5) 



where Nal„ are the observed values of the activity index, 
Nal„,^/„e are the activity values obtained from the fitted sine 
model, and (T/(Na l)„ the internal errors of each individual Na l„. 
Smaller p-values will indicate that the data was better fitted by 
a sinusoid. Note that these fits were carried out without taking 
in consideration whether the star passed the variability F-tests in 



Sect.[3]or not. Also, it is only possible to calculate the probabil- 
ities for the stars which have five or more data points, otherwise 
{N - A) < 0. Eight stars had their long-term activity success- 
fully fitted by a sinusoidal signal, with Pix^) < 0.1. These stars 
are Gil, G1382, G1433, G1581, G1667C, G1680, G1832, and 
HIP85647. However, Gl 1, G1382, and G1680 have not passed 
the long-term variability F-tests and therefore they will be dis- 
carded from discussion. The others will be discussed individu- 
ally in Sect.[6l 



5. Radial-velocity, activity and CCF parameters 

Figure [T] shows the time- series of radial velocity, activity, and 
the CCF parameters' BIS, FWHM and contrast for the 27 stars 
that passed the selection criteria (before the variability F-tests). 
Small points are nightly averaged data used to calculate the bins, 
points with errorbars are the binned data (the errors are explained 
in Sect. [2]). The peak-to-peak variation (A) and rms (cr) are also 
shown. The RV time- series have the known planetary compan- 
ions' signal subtracted. 

HIP38594 showed a trend in RV after removal of secular ac- 
celeration, with a slope of -50.83 ± 0.51 ms"^ yr"^. This cor- 
responds to a variation of A(y^) = 307.5 ms"^ in our timespan 
and is probably due to a stellar companion. A keplerian fit gives 
an orbit for a companion with a minimum period (only a linear 
trend is observed in the data) of 4941 ±1516 days. Since we are 
only concerned with low-amplitude and long-term variations, we 
removed this trend. 

Other stars appear to have linear trends in RV, such as Gl 205 
(with A(K) = 9 ms-i) and G1680 (A(K) = 15 ms'^), but since 
these trends have lower amplitude variations we decided to keep 
them for the rest of the study (see Sect. 15.21 for more informa- 
tion). 



5.1. Activity and radial-velocity scatter 

Table [2] shows the radial- velocity rms after the subtraction of 
the orbits together with the mean and rms of the Na I index for 
the binned data of the full sample. The scatter in radial velocity, 
cTeiVr), varies between 0.29 for Gl 849 and 5.9 m s"^ for Gl 680. 
The most active star in the sample is Gl 551 (Proxima Centauri) 
with an average Na l index of 0.49, while the most inactive star 
is Gl 699 (Barnard's star) with an average Na i value of 0.059 
(but highest relative scatter in Nai with -14%). We should alert, 
however, that the average Nai index values depend on stellar 
color (see Fig. 6 of Paper I) and should not be used to compare 
activity between stars with diff'erent colors without having the 
index calibrated for the photospheric contribution and without 
being reported to the bolometric flux. 

The most active star in the sample, G1551, also have the 
highest absolute scatter in Nai index (and 12% relative scatter). 
But its RV dispersion is not particularly high, having a standard 
de viation of (Tg(VV) = 1-5 m s"^ . 

iMeunier et al.l (l2010l) simulated the efl'ects of spots, plages, 
and inhibition of convection in the solar radial velocity during a 
solar activity cycle when the Sun is seen edge-on and observed 
as a star (integrated flux over the whole disk). The authors found 
that the signal induced by the three efl'ects on the solar RV during 
the cycle had an rms of cr(Vr) = 2 AO m s"\ with a peak-to-peak 
variation of 10.6 m s"^ . 

The median RV rms of 1.5 ms~^ f or our sample is lower 
than the results of IMeunier et al. l (l20TQh and therefo re in agree- 
ment with the extrapolation towards M dwarfs of iLovis et al.l 



4 



J. Gomes da Silva et aL: Long-term magnetic activity of a sample of M-dwarf stars from the HARPS program 

Table 2: Statistics for RV and Nai, and probabilities P(F) of the variability F-tests for activity. All parameters were calculated for 
the data with 150 day bins. 



Star Ntins Sel. o-^CK) < cr,(K) > <Nai> cr,(Nai) <cr,(Nai)> P{F) 

[ms~^] [ms~^] 



Gil 


5 


Yes 


0.58 


0.71 


0.104 


0.0021 


0.0013 


0.18 


G1176 


7 


No 


2.7 


0.94 


0.188 


0.0074 


0.0044 


0.12 


G1205 


6 


No 


2.9 


0.82 


0.179 


0.0046 


0.0029 


0.16 


G1273 


6 


Yes 


1.5 


0.49 


0.075 


0.0037 


0.0016 


0.042 


G1382 


5 


No 


2.1 


2.2 


0.196 


0.0052 


0.0027 


0.12 


G1393 


5 


No 


1.6 


0.71 


0.146 


0.0037 


0.0023 


0.18 


G1433 


6 


Yes 


1.2 


0.71 


0.141 


0.0034 


0.0011 


0.011 


G1436 


8 


Yes 


0.43 


0.52 


0.100 


0.0042 


0.0012 


0.0017 


G1479 


4 


No 


2.0 


0.97 


0.186 


0.0040 


0.0029 


0.31 


G1526 


4 


No 


1.8 


0.69 


0.158 


0.0035 


0.0011 


0.042 


G1551 


6 


No 


1.5 


0.83 


0.493 


0.057 


0.044 


0.28 


G1581 


13 


Yes 


0.96 


0.51 


0.071 


0.0026 


0.0013 


0.012 


G1588 


5 


Yes 


1.3 


0.39 


0.131 


0.0059 


0.0010 


0.0025 


G1667C 


8 


Yes 


0.73 


0.61 


0.096 


0.0030 


0.0014 


0.031 


G1674 


5 


No 


2.3 


0.73 


0.138 


0.0037 


0.0029 


0.32 


G1680 


6 


No 


5.9 


0.93 


0.146 


0.0048 


0.0022 


0.059 


G1699 


5 


No 


1.5 


0.61 


0.059 


0.0082 


0.0027 


0.027 


G1832 


5 


Yes 


0.30 


0.66 


0.132 


0.0039 


0.0015 


0.047 


Gl 849 


6 


Yes 


0.29 


0.77 


0.124 


0.0032 


0.0013 


0.032 


Gl 876 


6 


No 


2.3 


1.8 


0.081 


0.0047 


0.0017 


0.020 


Gl 877 


5 


Yes 


2.7 


1.3 


0.106 


0.0074 


0.0020 


0.014 


G1887 


4 


No 


1.6 


1.1 


0.173 


0.0038 


0.0014 


0.062 


G1908 


9 


Yes 


1.2 


0.55 


0.131 


0.0031 


0.0014 


0.020 


HIP12961 


6 


No 


1.3 


0.89 


0.172 


0.0028 


0.0016 


0.12 


HIP 19394 


7 


No 


5.3 


2.0 


0.094 


0.0033 


0.0019 


0.096 


HIP38594 


4 


No 


2.6 


1.2 


0.169 


0.0058 


0.0020 


0.055 


HIP85647 


6 


Yes 


1.5 


1.3 


0.145 


0.0041 


0.0017 


0.039 



(I2OIII) results which found that later-type stars have lower am- 
plitude RV influenced by long-term activity than earlier-type 
stars. It is expected that the stronger contribution to RV from 
magnetic activity cyc les comes from the inhibition of convection 
(iMeunier et al Since cell convection structure is strongly 

depending on the efl'ective temperature of the stars, the impact of 
their inhibition should be also depending of spectral type. Thus, 
it is expected that activity cycles will have difl'erent influence on 
RV for difl'erent types of stars. Furthermore, some of these stars 
may have undiscovered low-mass companions that could reduce 
the RV rms even further. 

Figure [2] shows the rms variation of radial velocity against 
the rms of Na i index (top panel) and relative variation of Na i 
(lower panel). A relation between these parameters is clearly vis- 
ible: stars with higher long-term scatter of activity also present 
higher long-term scatter in radial velocity. The minimum long- 
term relative variation we could detect in the Nai index was 
around 2%, with highest scatter reaching around 14%. G1699 
appears to be an outlier, more evidently when we consider the 
relative activity values. Although this star has an average value 
of criVr) when compared to the rest of the sample, it has the 
highest relative activity variation, almost the two times the value 
of the star with the second highest relative activity variation. 

5.2. Correlation between RV and activity 

Table[3]shows the correlation coefficients of the relation between 
RV and Na i with the respective FAPs. The FAPs were calculated 
using bootstrapping of the nightly averaged data, followed by a 
re-binning and a subsequent determination of the correlation co- 
efficient for each of the 10000 permutations (see also Paper I). 
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Fig. 2: Top: Radial- velocity scatter versus Nai scatter. Bottom: 
Radial- velocity scatter versus relative Na i scatter. 
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Table 3: Pearson correlation coefficients with respective FAPs. 



Star 


Vr vs Na I 


Nai vs BIS 


Nai vsFWHM 


Na I vs Contrast 




P 


FAP 


P 


FAP 


P 


FAP 


P 


FAP 


G1273 


0.84 


0.022 


-0.64 


0.094 


0.60 


0.13 


-0.89 


0.0074 


G1433 


0.91 


0.0018 


-0.70 


0.11 


0.70 


0.031 


-0.53 


0.14 


G1436 


0.83 


0.0091 


0.28 


0.29 


0.82 


0.0074 


0.44 


0.15 


G1526 


0.86 


0.13 


0.27 


0.34 


0.90 


0.057 


-0.76 


0.13 


G1581 


0.24 


0.25 


0.11 


0.38 


0.60 


0.019 


-0.38 


0.13 


G1588 


0.82 


0.033 


0.03 


0.49 


-0.51 


0.24 


-0.40 


0.25 


G1667C 


0.45 


0.12 


-0.35 


0.26 


-0.27 


0.22 


0.37 


0.17 


G1699 


0.32 


0.31 


-0.12 


0.42 


0.55 


0.15 


-0.79 


0.070 


G1832 


-0.20 


0.36 


-0.41 


0.30 


0.57 


0.19 


0.08 


0.46 


G1849 


-0.37 


0.25 


-0.38 


0.26 


0.69 


0.097 


-0.29 


0.29 


G1876 


0.63 


0.091 


0.45 


0.15 


0.27 


0.30 


-0.15 


0.39 


G1877 


0.77 


0.11 


-0.74 


0.085 


0.93 


0.023 


-0.45 


0.21 


G1908 


0.85 


0.0038 


-0.61 


0.051 


0.82 


0.0075 


0.43 


0.093 


HIP85647 


-0.03 


0.48 


0.32 


0.30 


0.81 


0.039 


-0.87 


0.014 



Notes. Only the 14 stars that passed the activity variability tests in Sect. [3] are considered. Bold highlights indicate FAP values lower than 0.05 
(95% significance level). FAPs were calculated via bootstraping as in Paper I. 



The tendency for positive correlations is clear, with an average 
correlation coefficient of < p >= 0.49. Five stars show a signif- 
icant correlation coefficient (with < 5% false-alarm probability) 
between activity and radial velocity, which represents 36% of 
the sample. These stars are G1273, G1433, G1436, G1588, and 
G1908 and their coefficients range between 0.82 and 0.91. One 
more star, Gl 876, has a FAP value lower than 10%, with a cor- 
relation coefficient of 0.63. If a FAP lower than 10% is taken as 
the limit, then we have ~43% of the variable stars having their 
RV induced by long-term activity. 



Fig. [3] shows the slope of the correlation between radial ve- 
locity and Nai index as a function of (V - /) colour. Symbols 
are in three sizes dependent on the correlation coefficient values: 
large for p > 0.75, medium for 0.50 < p < 0.75, and small for 
p < 0.50. The data points have three colors based on the FAP 
values: black for FAP < 0.01, grey for 0.01 < FAP < 0.05, and 
white for FAP > 0.1. The errorbars are the errors on the slope. 
We can see in the figure that there is a clear tendency for positive 
slopes of the correlation RV-Nai. This implies that a positive 
change in activity will induce a positive change in the apparent 
stellar velocity. 



All the cases of significant correlation have positive coeffi- 
cients and the slopes of the correlation range from 85 to 338 
m s"^ Na . Contrary to what lLovis et al ](|201I|) found for FGK 
stars, we found no strong relationship between the slopes and 
stellar colour (a proxy of Teff ) for this sample of early-M dwarfs. 
Furthermore, we found no cases of significant anti-correlations 
as could be expected from the author's study. It can be that 
the author's correlation observed in their Fig. 18 (lower panel) 
reaches a plateau for lower eflfective temperatures or that the 
variations induced by long-term activity on the RV of M dwarfs 
are so low that they are more difficult to observe (some of the 
observed scatter reaches the same level as the instrument pre- 
cision). Another possibility could be that still undetected small 
planets might interfere with the radial- velocities and therefore 
ruin any hypothetical correlation. 



cd 600 - 



OS 

I 



o 
o 




2.2 2.4 
V -I [mag] 

Fig. 3: Slope of the correlation between velocity and activity 
against (V-/) color. Symbol sizes depend on the value of the cor- 
relation coefficient between RV and activity: large for p > 0.75, 
medium for 0.50 < p < 0.75, and small for p < 0.50. The data 
points have three colors based on the FAP values: black for FAP 
< 0.01, grey for 0.01 < FAP < 0.05, and white for FAP > 0.1. 
The errorbars are the errors on the slope. 



5.3. Correlation between activity and the CCF parameters 
5.3.1. Nai vs BIS 

In this work we found a marginal tendency for negative correla- 
tions between Nai and BIS, with < p >= -0.18. No cases of sig- 
nificative coefficients with FAP < 5%. Three cases of marginal 
(5% < FAP < 10%) anti-correlation, G1273 with p = -0.64, 
G1877 with p = -0.74, and G1908 with p = -0.61 (21% of 
the sample). Due to the lack of significant cases of correlation 
with activity, the line bisector seems thus not to be a very good 
long-term activity indicator for early-M dwarfs. 

This is diflTerent from the trend found between activity and 
BIS for early- K stars where the tw o were found to be positively 
correlated (e.g. lSantos et al.ll2010l). While th e BIS values found 
here were all negative, ISantos etaP (l2010l) measured positive 
BIS values in all stars. Although the absolute value of BIS would 
increase with activity, for the case of negative values found in 
this study it means that the increase was a negative increase, and 
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therefore anti-correlated with activity. We can speculate that the 
difference in the signal of the BIS values might be attributed to, 
either the bisectores of K and M dwarfs having inverse shapes, 
or the use of differen t cross-correlatio n masks to obtain the CCF 
line profiles (see e.g. lDall et aDl2006h . 

5.3.2. Nai vs FWHM 

We found a tendency for positive correlations between our activ- 
ity indicator and width of the CCF profile, with < p >= 0.53 (the 
strongest average correlation coefficient between parameters). 
Similarly, a positiv e trend was also found for early-K stars (e.g. 
ISantos et al.ll201Ql) . And therefore we consider that the quali- 
tative behaviour of FWHM with activity i s similar fo r differ - 
ent spectral types, an effect also shown in iLovis et aP (l201ll) . 
Six stars have strong correlations in the range 0.60-0.93 with 
FAP < 1% (which represent 43% of our 14-star sample). These 
are G1433, G1436, G1581, Gl 877, G1908, and HIP85647. Two 
more stars show marginal correlation, G1526 and Gl 849. If we 
count them, then 57% of our stars have long-term correlation be- 
tween activity and FWHM. This means that the use of FWHM 
should be useful to detect long-term activity-like variations in 
M-dwarfs and can be used in complement to other activity prox- 
ies. 



5.3.3. Nai vs contrast 



As was detected by lSantos et al.l (l2010h for earlier- type stars, we 
found a tendency for anti-correlation, < p >= -0.30 between 
the activity level and the contrast of the CCF. But only two cases 
of significant coefficients, Gl 273 with p = -0.89 and HIP85647 
with p = -0.87 were found. This represents only 14% of the 
sample, an indication that the depth of the CCF line is not an 
optimal measure of the activity level of M dwarfs (at least as 
measured by the Nai/Caii lines). Two other cases of marginal 
correlation, Gl 699 with p = -0.79 and Gl 908 with p = 0.43 are 
present. 

6. Individual cases 

6.1. Stars with significant RV-activity correlation 

Gl 273 This star shows the highest slope of the RV-activity rela- 
tion of the stars with significant correlation. Its slope has a value 
of 338 m s"^ Na with the RV having a peak-to-peak variation 
of A = 5 ms"^ and a rms of ~1.5 ms"^ Its RV-activity corre- 
lation coefficient is 0.84 with a FAP of 2% (Fig. (H). This star 
also presents a strong anti-correlation between Na i and contrast 
with p = -0.89 (FAP = 0.0074). A moderate anti-correlation is 
also observed between Nai and BIS (p = -0.64, FAP = 0.094). 
Although we found signs of correlation between activity and RV 
we were not able to properly fit a sinusoidal to the Na i time se- 
ries. 



Gl 433 In Paper I we found this star to present a maximum in ac- 
tivity, typical of cycle-type long-term activity variations. Figure 
\5\ (middle and bottom panels) shows that the RV time- series, 
after removal of the planetary signal (Delfosse et al. 2012, in 
prep.), follows a similar pattern as the Nai index. The correla- 
tion coefficient between the two is p = 0.91 with a FAP of 0.00 1 8 
(Fig. [21 top panel). The Nai index also shows correlations with 
CCF parameters of -0.70 and 0.70 for BIS and FWHM respec- 
tively. We can therefore confirm that the activity cycle maximum 
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Fig. 4: Radial- velocity relation with Na I index for Gl 273. Small 
dots without errorbars are nightly averaged and points with er- 
rorbars are averaged over 150 days. The dashed line is the best 
linear fit. The correlation coefficient p and the respective FAP 
are shown. 



is inducing a maximum in RV. The slope of the RV-Na i corre- 
lation is 328 ms"^ Nai"^ and we measured an overall variation 
in RV of 3.29 ms"^ (with cr,(y,) = 1.22 ms"^). The ampli- 
tude and period of the cycle in both Na i or RV scales cannot be 
inferred because we do not have an entire cycle period in our 
data. However, a minimum limit for the amplitude and period 
can be calculated by fitting a sinusoid to the activity time- series. 
We determined that the minimum activity cycle period is 1665 
days with a minimum amplitude of 0.004 in Nai (Fig. [51 middle 
panel). A similar minimum period was also found for the radial 
velocity signal with a value of 1758 days and an amplitude of 
1.45 ms"^ (Fig. [3 bottom panel). This star represents therefore 
a good example of an RV signal induced by an activity cycle. 



Gl 436 This star has the smallest slope of the RV-activity cor- 
relation with a value of 85 ms"^ Nar^ From Fig. [6] we can ob- 
serve a large scatter in the nightly averaged data. Bal lard et"aD 
(I2OIOI) reported on short-term noise in photometry that they at- 
tributed to stellar spots. Furthermore, Knutson et al. (2011) de- 
tected evidence of occulted spots during the transits of G1436b 
(iButler et al.|[2004l) . It is probable that short-term activity vari- 
ability is producing the large scatter in Na i and contributing to a 
reduction of the slope's value. Nevertheless, the correlation co- 
efficient of the RV-activity relation is p = 0.83 with a FAP of 
0.9%, providing clues for the influence of long-term activity on 
the observed RV of the star. A strong correlation coefficient be- 
tween Nai and FWHM of 0.82 (FAP = 0.7%) was also found. 
However, we tried to fit a sinusoidal signal to both the RV and 
Nai time series but obtained no significant results. The activity 
time series in Fig. [T] shows a decreasing trend in activity that can 
be due to a long-period magnetic cycle, but more data will be 
needed to firmly prove this. 



Gl 588 With a correlation coeflScient of 0.82 (FAP = 3.3%) this 
star also presents a strong long-term relationship between ve- 
locity and activity (Fig. [7]). This correlation has a slope of 181 
ms"^ Nai"^. The RV time series (Fig.[T]) shows a periodic-like 
signal with a peak-to-peak variation of 3.6 ms"^ similar to the 
Nai time series apart from a point at BJD ~ 2453560 d. This is 
a point which includes only four nights that might be influenced 
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Fig. 5: Top: Radial- velocity relation with Nai index for G1433. 
Middle: sinusoidal fit to the Nai index for G1433. Bottom: 
Sinusoidal fit to RV of G1433. Small dots without errorbars are 
nightly averaged and points with errorbars are averaged over 150 
days. The dashed lines represent the best linear (top) and sinu- 
soidal (middle, bottom) fits to the data. 



by short-term activity variability. It is therefore not possible to 
conclude about the periodicity of this stars' activity cycle. 



Gl 908 This is another good example of long-term correlation 
between RV and activity. The correlation coefficient has a value 
of 0.85 (FAP = 0.4%) and the trend can be observed both in the 
binned and in the nightly averaged data (Fig. (8]). The slope of the 
correlation is 322 m s"^ Nai"^ and we measured an RV peak-to- 
peak variation of 3.85 ms"^. A strong correlation between Nai 
and FWHM with p = 0.82 (FAP = 0.75%) was also detected. 
Although we have strong evidence for RV induced by long-term 
activity, we could not fit a sinusoidal signal to both time series 
with significant pix^). This star has probably a magnetic "cycle" 
with no determined period, i.e., a star with a long-term aperi- 
odic magnetic activity variabi lity (these types of va riations were 
classified as var (variable) by lBaliunas et al.|[T995l) . 
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Fig. 6: Radial- velocity relation with Na I index for Gl 436. Small 
dots without errorbars are nightly averaged and points with er- 
rorbars are averaged over 150 days. The dashed line is the best 
linear fit. The correlation coefficient p and the respective FAP 
are shown. 
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Fig. 7: Radial-velocity relation with Na I index for Gl 588. Small 
dots without errorbars are nightly averaged and points with er- 
rorbars are averaged over 150 days. The dashed line is the best 
linear fit. The correlation coefficient p and the respective FAP 
are shown. 



6.2. Other interesting cases 

Gl 1 76 This star have what appears to be a long-term periodic 
signal with P ~2043 d in RV (see time-series plot Fig.[T]). The 
same trend is not observed in the Na i index and it did not passed 
the variability tests. This stars has a confirmed planetary com- 
panion, Gl 176 b with P = 8.7 d, and a rotationally modulated 
activity signal with a period of P = 39 d (Forveille et al. 2009b. 
But none of these signals can explain the -2000 d variation. It 
can be due to a long-period yet undiscovered planet. 

Gl 581 This is the star which has the best probability of having 
the activity data fitted by a periodic sinusoidal signal instead of a 
linear trend. We obtained a signal with a period of P = 3.85 years 
(P(X^) = 0.01%) and our timespan for this star is long enough to 
cover more than one activity cycle period (Fig.O. Furthermore, 
Gl 581 also passed the variability tests. To date, four planets are 
known to orbit the star with two more planets that were dis- 
carded (iBonfils et al. 2005: Udry et al. 2007; Mayor et al. 2009t 
iForveille et al .112011^ . We removed the signal of the four con- 
firmed planets which resulted in a long-term RV rms of 0.96 
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Fig. 8: Radial- velocity relation with Na I index for Gl 908. Small 
dots without errorbars are nightly averaged and points with er- 
rorbars are averaged over 150 days. The dashed line is the best 
linear fit. The correlation coefficient p and the respective FAP 
are shown. 
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Fig. 9: Sinusoidal fit to the activity time series of Gl 581. 



m s ^ which is of the order of the instrument precision. No corre- 
lation was found between long-term velocity and activity but the 
Nai index is correlated with FWHM (p = 0.60, FAP = 1.9%). 
Because of the large number of detected planets, it is sensible 
to recognise that their orbits might not have been properly sub- 
tracted from the data. We will therefore not take any conclusions 
about the possibility of the activity cycle being influencing the 
RV signal of this star. 
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Fig. 10: Sinusoidal fit to the activity time series of Gl 667 C. 



0.150 
0.145 
0.140 

I— I 

^ 0.135 
0.130 
0.125 



0.120 



G1832 




• 


" period 


= 1726.78 d (4.73 yr) 


• 






• • - 

• : 












• 




= 0.0323 





52500 53000 53500 54000 54500 
BJD - 2400000 [days] 



55000 



55500 



Fig. 11: Sinusoidal fit to the activity time series of Gl 832. 
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Fig. 12: Sinusoidal fit to the activity time series of HIP85647. 



star is moderately correlated with both FWHM (p = 0.55, FAP 
= 0.15) and contrast (p = -0.79, FAP = 0.070). 



Gl 6670 One planet is known to orbit this star (Delfosse et al. 
2012, in prep.). It is also a member of a triple system and or- 
bits the A -h B binary system. This system introduces a trend in 
RV that was subtracted together with the planetary companion 
signal. The Na i activity index passed the variability F-test and 
a sinusoidal function with a period ofP = 3.18 years is well 
fitted to the activity time series (with P(x^) = 3.3%, Fig. [TOl). 
However, the RV signal is only marginally correlated with ac- 
tivity, having a correlation coefficient of p = 0.45 with a FAP of 
12%. No correlations between the other parameters are detected. 

Gl 699 iKiirster et al.l (l2003l) found an anti-correlation between 
RV and the Ha index for the Bernard's star with a coeffi- 
cien t of p = -0.50. Using a longer timespan for the same 
star JZechmeister et al.l (l2009l) detected a similar correlation with 
p = -0.42. We found only a marginal positive correlation with 
the Nai index (p = 0.32, FAP = 0.32). The Nai index of this 



GI832 This star passed the activity variability tests but we 
found no correlation wit h radial velocity a fter the removal of 
the planetary companion (IB alley et al .1120091) . None of the other 
parameters are correlated apart from a marginal correlation be- 
tween activity and FWHM of p = 0.57 but with a high FAP of 
19%. The activity time series is well fitted by a sinusoidal func- 
tion with P(x^) = 3.2% and a minimum period of 4.73 years 
(Fig. [II]). 

H I P 85647 After removal of the signal of the planetary compan- 
ion (iForveille et al.ll20l"Tbb of this star and a linear trend due to 
a stellar companion, we found no correlation between RV and 
activity. There are strong correlations between Na i and the CCF 
parameters FWHM and contrast with p = 0.81 (FAP = 3.9%) 
and p = -0.87 (FAP = 1.4%), respectively. We successfully fit- 
ted a sinusoidal signal to the activity time series and obtained a 
period of 2.74 years (with P(x^) = 3.4%, Fig.fT2l). 
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7. Conclusions 

We used a sample of 27 MO to M5.5 dwarfs to study the relation- 
ship between long-term activity, radial velocity, and parameters 
of the cross-correlation function given by the HARPS pipeline. 
As the indicator for activity level we used the Na i index as sug- 
gested in Paper I. We binned the data to 150-day bins to average 
out high-frequency noise and removed any RV signals induced 
by known stellar or planetary companions. 

A selection of stars having long-term activity variability was 
carried out by using F-tests. This resulted in a subsample of 14 
variable stars, which means that around half of the stars in our 
sample have significant long-term activity variations. 

The activity time series of some stars were fitted by sinu- 
soidal functions to infer the activity cycle's period or minimum 
period. Five stars could be statistically well fitted by sinusoids 
and the inferred periods varied between 2.8 and 4.7 years. We 
should note, however, that our data covered close to one period 
or less of the fitted signals and therefore the cyclic nature of 
the signals cannot be fully established. Even if these stars have 
cyclic activity variability the periods obtained by the sinusoidal 
fitting should be regarded as minimum periods since the times- 
pan is not long enough to cover two periods. 

We obtained long-term velocity rms in the range 0.30 to 5.9 
ms"^ and a relative scatter in Nai index for the 14-star variable 
subsample was in the range -2-14%. These two parameters, 
cr(Vr) and (T(Nal), appear to be correlated. The median RV rms 
of 1 .5 m s"^ that can possibly be induced b y activity is lower than 
the one obtained for the case of the Sun by lMeunier et al.l (1201 Ol) 
as predicted by extrap olation of the trend in Fig. 1 8 (lower panel) 
of iLovis et al.i(l201 ll) for late-type stars. 

We then searched for correlations between RV, activity and 
the CCF parameters BIS, FWHM, and contrast. The main gen- 
eral results we obtained can be summarised as: 

- We found a general tendency for positive correlations be- 
tween long-term activity and radial velocity variations. No 
relation between the slope of the correlation and stellar color 
was found. 

- Five out of 14 stars with long-term variability have signif- 
icant correlation between activity and radial velocity. This 
amounts to 36% of our subsample. These stars are G1273, 
G1433, G1436, G1588, and G1908. The maximum peak-to- 
peak RV variation we obtained for stars with significant cor- 
relation between long-term activity and RV was ~5 ms"^. 
Only G1433 had its activity well fitted by a sinusoidal sig- 
nal even though the signal spans less than a full period. We 
determined a minimum period of 4.6 years for the activity 
cycle of this star, and similar minimum period of 4.8 years 
was found for the RV signal. The other stars could have ape- 
riodic activity cycles but more data is needed to confirm this 
hypothesis. 

- Although we found a general tendency for BIS to be anti- 
correlated with activity, those correlations were not statisti- 
cally significant. Only 21% of the stars with long-term vari- 
ability showed marginal correlations between activity and 
BIS. While we found a general trend for anti-correlations 
between activity and BIS for M dwarfs, Santos et al. (2010) 
found positive correlations between the S mw activity index 
and BIS for early-K stars. This is a curious result that should 
be investigated further. 

- We found that the long-term activity level was in general 
well correlated with FWHM. 43% of the stars with long- 
term variability had their activity level significantly corre- 
lated with FWHM. This is a good indication that the width 



of the CCF profile can be used to follow long-term activity 
in M dwarfs, as a complement to activity indices. A similar 
re sult was found prev iously for the case of early-K dwarfs 
bv lSantos et al.l (l2010h . 
- When compared with activity, the contrast showed a ten- 
dency for negative correlations with only 14% of significant 
cases. Th e tendency for anti-correlations was also found by 
ISantos et al . ( 2010) for early-K dwarfs, but since the fraction 
of significant coefl&cients is very low, this parameter does not 
seem to be a strong tool to diagnose long-term activity vari- 
ability. 

We found some cases of radial velocity variations induced 
by long-term activity. However, this does not seem to be a gen- 
eral trend. This can be due to several reasons including the bad 
frequency sampling of observations (since this is a program ded- 
icated to the search of planets, most of the observations were 
focused on short-term RV variability to find short-period plan- 
ets or with bad high-frequency sampling when the aim was to 
detect long-period companions), short timespan of the data (al- 
though we found some activity cycle periods to be of the 3 -year 
order some could have periods longer than the ~ 6-year span of 
our observations), and the fact that some of the stars might prob- 
ably have low-mass and/or long-period planets not yet detected 
(and thus their RV signal could be adding confusion to our ac- 
tivity study). A longer timespan of observations will certainly 
contribute to a better understanding of how long-term activity 
cycles influence the detected velocity of these late-type stars. 

In light of these results we must advise planet hunters to 
carefully check for long-term activity variations when faced with 
long-timespans of RV data since activity cycles could be adding 
noise or even mimicking the low-amplitude (<5 ms"^) signals 
of planetary companions. 
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Fig. 1: Time-series of radial- velocity, Nai index, and BIS, FWHIVI, and contrast of the CCF line profile. 
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Fig. 1 : Continued. 
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Fig. 1 : Continued. 
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Fig. 1 : Continued. 
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